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Abstract 

The first-order eikonal approximation is frequently adopted in interpreting the re- 
sults of A(e, e'p) measurements. Glauber calculations, for example, typically adopt 
the first-order eikonal approximation. We present an extension of the relativis- 
tic eikonal approach to A{e, e'p) which accounts for second-order eikonal correc- 
tions. The numerical calculations are performed within the relativistic optical model 
eikonal approximation. The nuclear transparency results indicate that the effect of 
the second-order eikonal corrections is rather modest, even at Q 2 ~ 0.2 (GeV/c) 2 . 
The same applies to polarization observables, left-right asymmetries, and differen- 
tial cross sections at low missing momenta. At high missing momenta, however, the 
second-order eikonal corrections are significant and bring the calculations in closer 
agreement with the data and/or the exact results from models adopting partial- wave 
expansions. 
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1 Introduction 



The eikonal approximation [1-3] has a long history of successful results in 
describing scattering processes like nucleon-nucieus scattering, heavy-ion col- 
lisions, and electroinduced nucleon-knockout reactions. The latter class of re- 
actions, usually denoted as A(e, e'p), provide access to a wide range of nuclear 
phenomena like short- and long-range correlations, relativistic effects, the tran- 
sition from hadronic to partonic degrees of freedom, and medium modifications 
of nucleon properties. The interpretation of A(e, e'p) data heavily relies on an 
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accurate description of the effect of the final-state interactions (FSI), i.e., the 
interactions of the ejected proton with the residual nucleus such as rescatter- 
ing and/or absorption. The eikonal approximation has been widely used to 
treat these distortions, either in combination with optical potentials [4-7], or 
with Glauber theory, its multiple-scattering extension [8-15]. 

The eikonal scattering wave functions are derived by linearizing the continuum 
wave equation for the ejected proton. Hence, the solution is only valid to first 
order in 1/k, with k the proton's momentum, and the eikonal approximation 
is suited for the description of reactions at sufficiently high energies. To extend 
the applicability to lower energies, Wallace [16] has developed systematic cor- 
rections to the eikonal scattering amplitude. Several authors have investigated 
the effect of higher-order eikonal corrections in elastic nuclear scattering by 
protons, antiprotons, and a particles [17,18], heavy-ion collisions [19-22], and 
inclusive electron-nucleus scattering [23]. The aim of this Letter is to deter- 
mine the influence of higher-order eikonal corrections on A(e, e'p) observables. 
To this purpose, we extend the relativistic optical model eikonal approxima- 
tion (ROMEA) A(e, e'p) framework of Ref. [7]. Our formalism builds upon the 
work of Baker [24], where an eikonal approximation for potential scattering 
was derived to second order in 1/k. Here, this work is extended to include the 
effect of the spin-orbit potential. 

The outline of this Letter is as follows. In Section [21 the second-order eikonal 
correction to the ROMEA model is derived. Section [3] presents the results of 
the A(e, e'p) numerical calculations. We look into how the second-order eikonal 
correction affects more inclusive quantities like the nuclear transparency, as 
well as truly exclusive observables such as the induced normal polarization 
P n , the left-right asymmetry Alt, and the differential cross section. Finally, 
in Section HI we state our conclusions. 



2 Formalism 

For the description of the A(e, e'p) reaction, we adopt the impulse approxi- 
mation (IA) and the independent-nucleon picture. Within this approach, the 
basic quantity to be computed is the transition matrix element [25] 

(J") = / dr¥^ s (r) >(f) e^0 Ql (f) . (1) 
Here, <b a , and ^ are the relativistic bound-state and scattering wave func- 

1 k,m s 

tions, with a\ the quantum numbers of the struck proton and k and m s the 
momentum and spin of the ejected proton. The relativistic bound-state wave 
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function is obtained in the Hartree approximation to the a — u model [26] 
with the Wl parametrization for the different field strengths [27]. The scat- 
tering wave function W~ ^ appears with incoming boundary conditions and is 

related to by time reversal. Furthermore, J M is the relativistic one-body 

k,m a 

current operator. Throughout this Letter, we use the Coulomb gauge and the 
CC2 form of > [28]. 

We now turn our attention to the determination of the scattering wave func- 
tion ^jj*^ • We start by considering the Dirac equation for a proton with 

relativistic energy E = \Jk 2 + Mjy and spin state \m s ^ subject to Lorentz 
scalar and vector potentials V s (r) and V v (r). The Dirac equation for the four- 
component spinor V^j^ (r) is converted to a Schrodinger-like equation for the 



upper component (r) [7,29] 
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+ V c {r) + V so (r) (o-L-ir-p) 
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The central V c (r) and spin-orbit V so (r) potentials are defined in terms of the 
scalar and vector ones, V s (r) and V v (r). The lower component w 



related to the upper one through 



(+) 
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r is 
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E + M N + V s (r)-V v (r) 



- - (+) r-\ 
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When solving Eq. ([2]) in the eikonal approximation, a standard procedure is 
to replace the momentum operator p by the asymptotic momentum k in the 
spin-orbit (V so (r) o ■ L) and Darwin (V so (r) (—if ■ p)) terms, as well as in the 
lower component ([3]). In literature, this is usually referred to as the effective 
momentum approximation (EMA) [30]. For the upper component, one puts 
forward a solution of the form 



(4) 



i.e., a plane wave modulated by an eikonal factor r](r). Here, iV is a normal- 
ization factor. 

In the ROMEA approach [7,29], which adopts the first-order eikonal approx- 
imation, Eq. (j2J) is linearized in p leading to a solution for the eikonal factor 
of the form 



77 ROMEA (f) = V ROMEA (b, z) = exp I -i 



M 



N 



dz'V opt (b,z') 



(5) 
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where r = (b, z) , the z axis lies along the momentum k of the proton, and 
V opt (b,z) = V c (b,z) + V so (b,z) (a ■ b x k — ikz). Despite the fact that it is 
written as an exponential phase, the solution (jHJ) is only valid up to first order 
in V opt /k. 

In what follows, we will derive an expression for the eikonal factor r/(r) that 
is valid up to order V opt /k 2 . The momentum dependence in the spin-orbit 
and Darwin terms makes that these terms are retained up to order V so /k, 
while central terms are included up to order V c /k 2 . Note that the expansion is 
not expressed in terms of the Lorentz scalar and vector potentials V s and V v . 
Looking for a solution of the form (Tj0) for the Schrodinger-like equation (j2J), 
Baker arrived at the following equation for the eikonal factor (see Eq. (14) of 
Ref. [24]): 



T,(b,z) = l-i-jZ J dz'V opt (b,z') V (b,z') + -^ Vopt (b,z)rj(b,z) 

— oo 

— oo x ' 

Note that, apart from dropping contributions of order V opt /k 3 and higher, no 
additional assumptions were made when deriving Eq. (jSj). In Ref. [24], Eq. (jSJ) 
was subsequently solved for spherically symmetric potentials. The spin-orbit 
and Darwin terms, however, break the spherical symmetry and a novel method 
to solve Eq. is needed. To that purpose, we assume that the derivative of 
the function r\ is of higher order in 1/k than rj itself (as is true for the ROMEA 
solution (j2D). This allows us to drop the drj/db contribution in the last term 
of Eq. (jSJ), as it is of order V opt /k 3 or higher: 



|f / dz ' (z - z') (\ + 1) | (y opt (6, z ') r,(b, z ') 
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db 



n(b,z'). 



(7) 



Spherical symmetry implies that z' dV c (b, z')/db = bdV c (b,z')/dz'. Hence, 
the z' dV c /db term in Eq. ([7D can be written as 
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In the first step, we made use of the fact that the derivative dr]/dz' is of higher 
order to turn the integrand into an exact differential. A similar reasoning, 
followed by integration by parts, leads to 



2k 2 \b db. 
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for the Darwin term of Eq. (j7j). Inserting the expressions of Eqs. 
Eq. (jnD adopts the form 
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We look for a solution of the form 



M 



v(b,z) = f(b,z) exp I -i ^ / dz' V opt (b, z') f(b, z') 



f(b,z) exp (iS(b, 



(11) 
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which should reduce to the ROMEA result of Eq. (JSJ) when terms of higher 
order than V opt /k are neglected. Accordingly, the function f(b, z) should be of 
the form / = 1 + 0(V opt /k 2 ). Substituting ffTTj) into Eq. (jTDJ) and multiplying 

by e~ lS( - b ' zS) yields 



f(b,z) 
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In deriving this equation, we set e* 5 '- 6 ' 2 '*' e * 5 ( b ' z ) equal to 1, since higher-order 
terms are neglected. The difficulty in solving for f(b, z) is that Eq. ( Tl2l) is an 
integral equation. An expression for f(b,z) can, however, be readily obtained 
by adding (1 — /) terms, which introduce only higher-order terms, to the right- 
hand side of Eq. ffl2l . This is permitted since we seek for a solution up to order 
V op t/k 2 . With this manipulation, the function / becomes 



2fc 2 V 96 / v ' y 2fc 2 6 V 96 J 96 



-oo 



+ 2^ Ko(6, z)(a-bxk- ikz) 

— oo 

/ dz' f2 + 6|-l ^o(M) • (13) 



2A: 7 \ 96, 

—00 v ' 

The eikonal factor of Eq. ffTTl) with f(b,z) given by (|13[) . is a solution of the 
integral equation (jBJ) to order K P tA 2 and reduces to the ROMEA result (jSJ) 
when truncated at order V^ p t/&. Furthermore, it can be easily verified that the 
derivative of r\ is of higher order in V opt /k than r\ itself. Henceforth, calculations 
performed with the eikonal factor of Eqs. (TTTT) and (|T3l) . are dubbed as the 
second-order relativistic optical model eikonal approximation (SOROMEA). 
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3 Results 



One way to quantify the overall effect of FSI in A(e, e'p) processes is via 
the nuclear transparency. The measurements are commonly performed under 
quasielastic conditions [31-36]. We obtain the theoretical transparencies by 
adopting similar expressions and cuts as in the experiments. Hence, the nuclear 
transparency is defined as [37] 

T = E a S& Pm dp m S a {p m , E m , k) 

C A J2a lA 3 p m dpmSpyj lA (p m , E m ) 



Here, S a is the reduced cross section for knockout from the shell a 

dh ° a ( e e'v) 

S a (p m , E m , k) = dn ^' ? , (15) 

Ka ep 



where p m and E m are the missing momentum and energy, K is a kinemati- 
cal factor and a ep is the off-shell electron-proton cross section. 5"p WIA is the 
reduced cross section within the plane-wave impulse approximation (PWIA) 
in the nonrelativistic limit. Further, J2 a extends over all occupied shells a in 
the target nucleus. The phase-space volume in the missing momentum A 3 p m 
is defined by the cut \p m \ < 300 MeV/c. The A-dependent factor c A corrects 
in a phenomenological way for the effect of short-range correlations. We intro- 
duce the ca in the denominator of Eq. (fUj) because the data have undergone 
a rescaling with c A = 0.9 ( 12 C) and 0.82 ( 56 Fe). 

Transparencies have been computed for the nuclei 12 C and 56 Fe at planar 
and constant (q, u) kinematics compatible with the phase space covered in 
the experiments. For the optical potential, the EDAD1 parametrization of 
Ref. [38] was used. 

In Fig. [1] the ROMEA and SOROMEA results are displayed as a function 
of the four-momentum transfer Q 2 and compared to the data. Not surpris- 
ingly, at high Q 2 , the ROMEA and SOROMEA predictions practically coin- 
cide and the role of the second-order eikonal effects grows with decreasing Q 2 . 
At Q 2 = 1.7 (GeV/c) 2 , the ROMEA and SOROMEA transparencies agree to 
within 1%; while at Q 2 = 0.3 (GeV/c) 2 , the difference has risen to 3% for 
56 Fe and 5% for 12 C. The enhancement of the nuclear transparency due to 
the second-order eikonal corrections is modest, even for values of the four- 
momentum transfer as low as Q 2 = 0.2 (GeV/c) 2 . Both the ROMEA and 
the SOROMEA predictions tend to slightly underestimate the measurements. 
The second-order corrections move the predictions somewhat closer to the 
Q 2 = 0.34 (GeV/c) 2 data point. 
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As the nuclear transparency involves integrations over missing momenta and 
energies, it may hide subtleties in the theoretical treatment of the FSI mech- 
anisms. Next, we focus on highly exclusive A(e, e'p) quantities and quantify 
the role of second-order eikonal effects. 

An observable that is particularly well suited to study FSI effects is the induced 
normal polarization 

p = d 5 a(a n =l)-d 5 <T K=I) 
d 5 a (a n =T) + d 5 a (a n = j) 



where a n denotes the spin orientation of the ejectile in the direction orthogonal 
to the reaction plane. Indeed, in the one-photon exchange approximation, P n 
vanishes in the absence of FSI. 

Fig. [2] shows the missing momentum dependence of the induced normal polar- 
ization for the kinematics of Ref. [39], corresponding with Q 2 « 0.5 (GeV/c) 2 . 
The calculations are performed with the energy-dependent A-independent 
(EDAI) potential of Ref. [38]. The ROMEA results are in line with the rela- 
tivistic distorted- wave impulse approximation (RDWIA) calculations of Ref. [40] 
The RDWIA framework was implemented by the Madrid- Sevilla group [41] 
and relies on a partial-wave expansion of the exact scattering wave function. 
It is similar to the (SO)ROMEA approach in that both models compute the 
effect of the FSI with the aid of proton-nucleus optical potentials. Further, 
the overall agreement with the data is excellent. The second-order eikonal 
corrections are most pronounced for the lsi/2 level. For missing momenta 
p m > 125 MeV/c, they reduce the magnitude of the P n for the lsi/ 2 state 
by roughly 20%, thereby resulting in a marginally better agreement with the 
highest p m data point. For lp^/2 knockout, on the other hand, the effect of the 
second-order eikonal corrections is smaller than 5%. 

The inclusion of the second-order eikonal effects is particularly visible at high 
missing momentum, a region where also other mechanisms become impor- 
tant. The qualitative behavior of the meson-exchange and A-isobar currents, 
for instance, is alike [42]. At low missing momenta (p m < 200 MeV/c), the in- 
duced normal polarization P n is relatively insensitive to the two-body currents; 
whereas at higher missing momenta, sizable contributions from the meson- 
exchange and isobar currents are predicted. The influence of the meson and 
isobar degrees of freedom is also stronger for knockout from the lsiy 2 shell 
than for lp^/2 knockout. 

In Fig. [21 also calculations neglecting the spin-orbit part V so (b, z) 3 -b x k are 
shown. They illustrate that the spin-orbit distortion is the largest source of 
P n . Hence, a correct inclusion of this term is essential. Moreover, P n proves to 
be rather sensitive to the choice of optical potential [40]. 
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Another A(e, e'p) observable which has been the subject of many investigations 
is the left-right asymmetry 

_ d 5 a (0 = 0°) -^a (0 = 180°) 
LT d 5 a (0 = 0°) + d 5 a (0 = 180°) ' 1 ' 



The subscript LT indicates that this quantity is closely related to the longitudinal- 
transverse response function. 

Fig. [3] presents the Alt predictions for the removal of lp-shell protons in 
16 in the kinematics of Refs. [43,44]. The FSI shift the dip in Alt, which 
is located at p m « 400 MeV/c in the relativistic PWIA (RPWIA), to lower 
values of the missing momentum. This shift is essential to describe the data at 
p m pa 350 MeV/c. The exact p m location and height of the ripple, however, are 
affected by many ingredients of the calculations, such as the current operator, 
bound-state wave function, and parametrization of the optical potential [44]. 
As can be inferred from Fig. [3j the second-order eikonal corrections affect the 
height, but not the position of the ripple. 

We also show the results of our SOROMEA calculations within the so-called 
noSV approximation. In this approximation, the dynamical enhancement of 
the lower component of the scattering wave ([3]) due to the V s (r) — V v (r) term 
is omitted. As such, the SOROMEA-noSV calculations make the same set of 
assumptions as the EMAf-noSV predictions by the Madrid-Sevilla group. The 
EMAf-noSV approach is an RDWIA calculation which adopts the EMA in 
combination with the noSV approximation. The second-order eikonal correc- 
tions clearly increase the height of the oscillation in Alt and brings the eikonal 
noSV calculations in excellent agreement with the corresponding partial-wave 
prediction EMAf-noSV. Finally, the comparison between the SOROMEA and 
the SOROMEA-noSV calculations demonstrates that the dynamical enhance- 
ment plays a significant role in the description of the Alt data. 

In Fig. HI 16 0(e, e'p) cross-section results are displayed for the kinematics 
of Fig. [31 The spectroscopic factors, which normalize the calculations to the 
data, were determined by performing a \ 2 fit to the data and are summa- 
rized in Table [TJ The RDWIA spectroscopic factors are 5-10% higher than 
the (SO)ROMEA ones. The second-order eikonal corrections hardly affect the 
values of the extracted spectroscopic factors. Both our (SO)ROMEA calcula- 
tions and the RDWIA predictions of the Madrid-Sevilla group do a very good 
job of representing the data over the entire p m range. For missing momenta 
\p m \ < 250 MeV/c, the (SO)ROMEA and RDWIA results are in excellent 
agreement. The impact of the second-order eikonal corrections on the com- 
puted differential cross sections is almost negligible for p m below the Fermi 
momentum, but can be as large as 30% at high p m . The inclusion of the 
second-order effects improves the agreement with the RDWIA calculations at 
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these high missing momenta. Results for the effective response functions Rl, 
Rt, Rlt-, and Rtt are not shown, but the effect of the second-order eikonal 
corrections is similar to the effect on the differential cross section. 





RPWIA 


ROMEA 


SOROMEA 


RDWIA 


1P3/2 


0.55 


0.84 


0.83 


0.92 


lPl/2 


0.47 


0.75 


0.74 


0.78 



Table 1 

The spectroscopic factors for the 16 0(e, e'p) reaction of Ref. [43], as obtained with 
a x 2 procedure. 



4 Conclusions 

We have developed a formalism to account for second-order corrections in the 
eikonal approximation. Our model is relativistic and includes both the central 
and spin-orbit parts of the optical potentials. The formalism has been applied 
to A(e, e'p) processes. Our numerical calculations show that the effect of the 
second-order eikonal corrections on A(e, e'p) observables is rather limited for 
Q 2 > 0.2 (GeV/c) 2 . The nuclear transparency calculations confirm the ex- 
pected energy dependence of the eikonal corrections: the effect decreases with 
increasing Q 2 . Concerning the p m dependence of the A(e, e'p) observables, the 
effect of the second-order eikonal corrections is minor except at high missing 
momenta. In this high-p m region, the eikonal corrections affect the observables 
up to an order of 30%, thereby bringing the calculations closer to the data 
and/or the RDWIA calculations. The robustness of the first-order eikonal ap- 
proximation, which emerges from this study, can be invoked to explain the 
success of the Glauber approach to A(e, e'p) down to relatively low kinetic 
energies of 200 MeV. 
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Fig. 1. Nuclear transparencies versus Q 2 for A(e,e'p) reactions in quasielastic kine- 
matics. The SOROMEA (dashed lines) are compared to the ROMEA (solid lines) 
results. The EDAD1 potential [38] has been employed in both formalisms. Data 
points are from Refs. [31] (open squares), [32,33] (open triangles), [34,35] (solid 
triangles), and [36] (open diamonds). 
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Fig. 2. Induced normal polarization P n for proton knockout from the 1^3/2 (upper 
panel) and lsi/2 (lower panel) shell in the 12 C(e, e'p) reaction. The kinematics is de- 
termined by beam energy e = 579 MeV, momentum transfer q = 760 MeV/c, energy 
transfer u = 292 MeV, and azimuthal angle (ft = 180°. The solid (dashed) curves 
represent ROMEA (SOROMEA) calculations. The dot-dashed (dotted) curves re- 
fer to predictions obtained within the ROMEA (SOROMEA) frameworks, with the 
spin-orbit term V so (b, z) a ■ b x k turned off. The data are from Ref. [39]. 
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Fig. 3. The left-right asymmetry Alt for the 16 0(e,e'p) experiment of [43]. 
The kinematics was e = 2.442 GeV, q = 1 GeV/c, and u = 445 MeV (i.e., 
Q 2 = 0.8 (GeV/c) 2 ). The red solid (green dashed) lines show the results of the 
ROMEA (SOROMEA) calculations. The SOROMEA-noSV (orange long-dotted 
curves) calculations differ from the SOROMEA calculations in that the dynamical 
enhancement of the lower component of the scattering wave function is neglected. 
The cyan short-dotted curves present the results from an RDWIA calculation where 
the spinor distortions in the scattered wave are neglected. All calculations use the 
EDAI version for the optical potentials [38]. The black short-dotted curves represent 
the RPWIA results. The data points are from Ref. [43]. 
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Fig. 4. 16 0(e,e'p) cross sections compared to ROMEA, SOROMEA, RDWIA, and 
RPWIA calculations for the constant (q, oj) kinematics of Fig. [3l The calculations 
use the optical potential EDAI [38]. The data are from Ref. [43] and the RDWIA 
results from Ref. [44]. The following convention is adopted: positive (negative) p m 
corresponds to 4> = 180° ((f) = 0°). 
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